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Abstract 

We discuss a new approach to uhrarelativistic heavy-ion colhsions based on clas- 
sical hadrodynamics for extended nucleons, corresponding to nucleons of finite size 
interacting with massive meson fields. This new theory provides a natural covariant 
microscopic approach that includes automatically spacetime nonlocality and retar- 
dation, nonequilibrium phenomena, interactions among all nucleons and particle 
production. In the current version of our theory, we consider extended unexcited 
nucleons interacting with massive neutral scalar (o") and neutral vector (uj) meson 
fields. The resulting classical relativistic many-body equations of motion are solved 
numerically without further approximation for soft nucleon-nucleon collisions at 
Plab = 14.6, 30, 60, 100 and 200 GeV/c to yield the transverse momentum imparted 
to the nucleons. For the future development of the theory, the isovector pseu- 
doscalar (vr"'', vr", vr'^), isovector scalar ((5+, 5~ , 5^), isovector vector p~ , p^) 
and neutral pseudoscalar (r/) meson fields that are known to be important from 
nucleon-nucleon scattering experiments should be incorporated. In addition, the 
effects of quantum uncertainty on the equations of motion should be included by 
use of techniques analogous to those used by Moniz and Sharp for nonrelativistic 
quantum electrodynamics. 

1. Introduction 

Many of you are involved in the search for the quark-gluon plasma — a predicted new 
phase of nuclear matter where quarks roam almost freely throughout the medium instead 
of being confined to individual nucleons. Experimental identification of the quark-gluon 
plasma will require accurate predictions for ultrarelativistic heavy-ion collisions on the 
basis of conventional nuclear physics, with the known hadronic degrees of freedom properly 
taken into account. Significant deviations between these predictions and experimental 
data would then signal the onset of new phenomena such as a quark-gluon plasma. 

With this goal in mind, we have developed at Los Alamos an entirely new microscopic 
many-body approach to ultrarelativistic heavy-ion collisions based on classical hadrody- 
namics for extended nucleons, corresponding to nucleons of finite size interacting with 
massive meson fields. This approach, which satisfies a priori the physical conditions that 
exist at ultrarelativistic energies, is manifestly Lorentz-covariant and allows for nonequi- 
librium phenomena, interactions among all nucleons and particle production. 

The physical input underlying this approach consists of Lorentz invariance (which 
includes energy and momentum conservation), nucleons of finite size interacting with 



massive meson fields and the classical approximation applied in domains where it should 
be reasonably valid. This starting point builds upon the traditional hadronic descrip- 
tion of nuclear processes that has in the past been so successful in a wide variety of 
situations. By treating the underlying quarks and gluons implicitly in terms of nucleons 
and mesons, we are able to solve the resulting classical relativistic many-body equations 
of motion numerically without further approximation. In particular, we do not need to 
make either a mean-field approximation, a perturbative expansion in coupling strength 
or a superposition of two-body collisions. 

The motivation for the classical approximation in this new theory is that at bombard- 
ing energies of many GeV per nucleon, the de Broglie wavelength of projectile nucleons 
is extremely small compared to all other length scales in the problem. In addition, the 
Compton wavelength of the nucleon is small compared to its radius, so that effects due to 
the intrinsic size of the nucleon dominate those due to quantum uncertainty. Finally, the 
angular momentum is typically several hundred h and the radiated energy corresponds to 
several meson masses. The classical approximation for nucleon trajectories should there- 
fore be valid, provided that the effect of the finite nucleon size on the equations of motion 
is taken into account. 

We describe in sect. 2 the present version of our theory, which includes the neutral 
scalar (a) and neutral vector (u) meson fields. This permits a qualitative discussion of 
such physically relevant points as the effect of the finite nucleon size on the equations of 
motion, an inherent spacetime nonlocality that may be responsible for significant collective 
effects and particle production through massive bremsstrahlung. The a and u mesons that 
are produced through this mechanism will subsequently decay primarily into pions with 
some photons also. The resulting classical relativistic equations of motion are solved in 
sect. 3 for soft nucleon-nucleon collisions at piab = 14.6, 30, 60, 100 and 200 GeV/c to yield 
the transverse momentum imparted to the nucleons. Section 4 discusses future directions 
for the svstematic development of the theory. Further details are given in a series of 
papers,0"0 although not all of the equations appearing in some of the earlier publications 
are in their final form. 



2. Equations of motion 

Our action for extended unexcited nucleons interacting with massive scalar and 
vector meson fields is 
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where Mq is the bare nucleon mass and qi is the four-position of the ith nucleon, whose 
trajectory is given by qi = qi^ri). A dot represents the derivative with respect to Tj. In the 
action the four- velocities are not constrained so that g-^ = 1 and Tj is not yet identified 
as the proper time; it is only in the equations of motion, which are derived as a result of 
the variation of /, that this is true. We use the metric g^^^ = diag(l, —1, —1, —1), write 
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four- vectors as = q) = (g*, q^, q^ ,q^) and use units in which h = c = 1. The scalar 
potential is denoted by 0, the four- vector potential by V and the meson masses by rus^v 
The vector field strength tensor is 



the scalar source density is 



3{x) = 9sJ2j d'^i - 'ii)\f^ (3) 
and the vector source density is 
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where p is the four-dimensional mass density of the nucleon, the spatial oart of which 
we assume to be exponential in the nucleon's rest frame. The valuesiiryj that we are 
currently using for the six physical constants appearing in our theory are nucleon mass 
M = 938.91897 MeV, scalar (cr) meson mass rus = 550 MeV, vector (a;) meson mass 
TTT-v = 781.95 MeV, scalar interaction strength gj^ = 7.29, vector interaction strength 
(7v^ = 10.81 and nucleon r.m.s. radius Rrms = 0.862 fm. 

In ref. |^ we have derived exact equations of motion for the above action in two limits: 
(1) relativistic point nucleons and (2) nonrelativistic extended nucleons. We then gener- 
alize covariantly to obtain relativistic equations of motion for extended nucleons, which 
can be written as 

~ fs,i + fv,i + /s^ext,i + fv,ext,i ) (5) 

where M* is the effective mass, f^^ is the scalar self-force, f^^ is the vector self-force, 
fsext i is the scalar external force and ext i is the vector external force. These equations 
of motion are second-order, nonlinear, integrodifferential equations with four dimensions 
per particle. 



3. Transverse momentum for soft nucleon- nucleon collisions 

To solve our equations of motion (^) we use a fourth-order Adams- Moulton predictor- 
corrector algorithm with adaptive step sizes. The integrations over proper time are done 
with a special error-minimizing application of Lagrange's four-point (cubic) interpolation 
formulas. We present here results on the transverse momentum imparted in the soft 
collision of two nucleons at laboratory momentum joiab = 14.6, 30, 60, 100 and 200 GeV7c. 
At three of these momenta substantial experimental data exist for heavy-ion collisions0'El 
and at the remaining two momenta experimental data exist for proton-proton coUisions.El 
Results calculated for the scattering angle, radiated energy and rapidity are described in 
ref. I 

As shown in fig. 0, for a given incident momentum, the transverse momentum for soft 
reactions has a maximum value at a certain impact parameter and decreases to zero for 
both head-on and distant collisions. The maximum transverse momentum increases slowly 
with increasing incident momentum in this range, and the impact parameter at which 
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Fig. 1: Calculated dependence of the transverse momentum upon impact parameter for soft nucleon- 
nucleon collisions at five incident laboratory momenta. 



this maximum occurs decreases. For ultrarelativistic collisions this impact parameter is 
approximately the distance at which the transversely dominating static vector force for 
extended nucleonsi'S has its maximum. At the other extreme of low incident momentum, 
the opposing scalar and vector forces are of similar magnitude and give rise for small 
impact parameter to the more complicated behaviour of the double-dot-dashed curve in 

The qualitative behaviour of our results can be understood in terms of the nature 
of the external forces. The repulsive vector force scales as the Lorentz factor 7 in both 
the longitudinal and transverse directions, whereas the attractive scalar force scales as 
7^ in the longitudinal direction and as unity in the transverse direction. This implies 
that the vector force will dominate the transverse acceleration and the scalar force will 
dominate the longitudinal acceleration. For a given impact parameter the scattering 
angle and transverse momentum will be essentially proportional to the vector interaction 
strength g^"^, and the radiated energy will be essentially proportional to 7 times the scalar 
interaction strength g^^. 

As an initial test of the theory to describe the gross features of soft nucleon-nucleon 
collisions, we compare in fig. |^ experimental and theoretical values of the mean transverse 
momentum <pt> as a function of total centre-of-mass energy ^/s. Because the classical 
total cross section diverges, it is necessary when calculating such average quantities to 
select a particular cutoff in impact parameter, or equivalently, a value of the total cross 
section, for which purpose we have usedB fJtot = 40 mb. Our results calculated with the 
values of the constants listed in sect. 2 are given by the solid curve in fig. ^ and are seen 
to be both lower in magnitude and to increase more rapidly with ^/s than the experi- 
mental values.0lli However, the theoretical curve is calculated for all values of rapidity 
corresponding to a 40 mb total cross section, whereas the experimental results are peak 
values in the central region, where there is a short plateau in rapidity. Furthermore, the 
experimental values include high-px contributions arising from hard collisions. Because 
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Fig. 2: Comparison of experimental and theoretical values of the mean transverse momentum <pt> a-s a 
function of total centre-of-mass energy ^/s. 

of the rough proportionahty of <Pt> to the vector interaction strength g^'^, we show with 
the dotted curve in fig. || a simple estimate obtained by multiplying the solid curve by the 
ratio of g^"^ = 17.26 obtained by Bryan and ScottEj from an analysis of nucleon-nucleon 
scattering at laboratory kinetic energies between and 350 MeV to the value appearing 
in sect. 2. 



4. Future directions 

From nucleon-nucleon scattering experiments we know that several additional meson 
fields are important and must be included for a quantitative description^ 

• Isovector pseudoscalar (7r+, n^, 7r°) 

• Isovector scalar (5+, 6^ , 6^) 

• Isovector vector , p^) 

• Neutral pseudoscalar (77) 

The next step in the systematic development of the theory should be the inclusion of these 
additional meson fields. In addition, the effects of quantum uncertainty on the equations 
of motion should be studied and included if they are important. This should be possible 
by use of techniques analogous to those used by Moniz and Sharp for nonrelativistic 
quantum electrodynamics. El 

Once these two formal items are completed, the full version of the theory should be 
used to calculate ultrarelativistic heavy-ion reactions that are being studied experimen- 
tally. Systematic comparisons between the calculated and experimental results for such 
quantities as differential cross sections, transverse momentum distributions, particle mul- 
tiplicity distributions, rapidity distributions and particle-particle correlations will enable 
us to decide the usefulness of classical hadrodynamics for understanding these phenomena. 

This work was supported by the U. S. Department of Energy. 
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